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ABSTRACT 

Rtnr'* Mouiar onaration of thB DOE/NASA Mod-l wind tuiOlne Degan In Octobar 1979 about 10 nwrby 

navT^coiwlalned of noise from tnt' iMcnine. Oovelopment of the NASA-ueRC turbine touf^ 
n^SSlon 2oS SSn in May 1980 as part of an effort to undoratanci and reduce tne nolso generated by 
Mod-l *^Tone sSund^levels predicted aitn this code are in generally good agreement witn measured data 
taken in tne vicinity Mod*l wind turbine (lass than 2 rotor diameters). Comparison in tne far field 
iSlcates tnat propagation effects due to terrain and atmospheric conditions may be amplifying 
sSnS l2« s ny Sbout 6 db. Parametric analysis using the code has sno-n that the proowlnan 
Ss t^Mod-K rotor noise are (1) the velocity deficit in the wake of the support tower, (2) tne high 
rotor speed, and (3) off-optimum operation. 


the actual 


INTROOUCTION 


Since regular operation of the Mod»l wind 
turbine began in October 1979 about ten 
householos have complained of noise from the 
machine. The character of tne noise is 
described as an audible "tnump" at a repetition 
rate eoual to the blade passing freouency. In 
some instances, low freouency acoustic energy 
has resulted in complaints of vibrations within 
the homes. 

Since January 1980, efforts have been directed 
at identifying tne causes of the noise and 
methods of reducing it. One effort has been the 
development of computer prediction c^de5 and 
ootfiinlng experimental data to verify these 
codes. The objective of this effort is to 
unoerstand the noise generated by Mod-1 and to 
prevent it from being a problem on advanced 
machines. Development of a wind turbine sound 
prediction code began at the NASA-LaRC in May 
1980. This report presents (1) measured dots 
thet characterize ths noise problem at Mod-1, 

(2) the analytical method used in the NASA LeRC 
code, (3) validation of the code using experi- 
mental data from Mod-1, and (A) applications of 
the coiie to tne Mod-i situation. 

SUMMARY 

Mf'asureu data taken at the MOd-l site show that 
the im^julsive cnaracter of tne noise is composfib 
of hamonics of tne blade passing freouency. 

While some of these hai'monlcs exist below the 
audible t’reduetK^y range (about 20 to 20, OCX) Hz), 
many are aoove the nominal audible threshold of 
20 HZ. 

A computer program (WTliOUNU) has been developed 
for calculating the intensity and freouency 
chaiacterlstlcs of sounds generated by wind 
turbines in a non-uniform wind flow field. The 
results calculated with this code are generally 
in good agreement with Mod-l measured sound spec- 
tra. However, propagation effects due to terrain 
and atmo8[)»'eric conditions have complicated the 
amplitude correlation with Mod-l data in the far 
field. These effects have been estimated by the 
code to cause an amplification of 6 dB or more 
at a home that has registered complaints. 


The WTSCXJND code has been used to determine the 
source of the noise generation from moo-1 as 
well as identifying operating conditions 
associated with the hlgf^est noise levels. 

Mooellng of the wind flow field characteristics 
shows that the predominant source of noise from 
Mod-1 is the wind velocity deficit in the wake 
of the tower. Because the rotor plane is 
downwind of the tower, this deficit produces 
changes in the aerodynamic forces on the blades 
resulting in sound pressure variations in tne 
acoustic field. The level of the sound pressure 
variations is most strongly a function of rotor 
speed and wind speed. Reducing the rotor speed 
from 35 to 23 rpm is predicted to reduce souno 
levels by about 11 dB. The Increase in sound 
levels with wlndspeed is expected to be about 12 
(fi between cutin and rated windspeeds. vari- 
ation in sound level with power is not clear-cut 
since high sound levels can occur during both 
high and low power conditions. 

MOD-1 )<EASURED DATA 

In January 1980, the Solar Energy Research 
Institute (SERI) and the General Electric 
Company under DOE funding began visiting the 
Mod-l site to investigate the noise problem. 
During thoie visits, tape recordings of tne 
sound pressures were made botn near the Mod-l 
and at residences that had complained of noise. 

A good example of the recorded time history of 
the pressure variation near the Mod-l recorded 
by SERI is shown in Elgure 1. The Impulsive 
nature is clearly apparent e$ eacn blade passes 
behind the tower structure. The passage of the 
blade throu^ the wake of each leg of the tower 
is seen as the two upper peaks approximately 
0.075 seconds apart. 

To better understand the character of the 
pressure Impulse, spectrum analyses were 
performed on the measured data, One sample 
spectrum of the sound near the Mod-l is showt* In 
Figure 2. This was recorded v^'en the machine 
was generating a heavy "thumping" sound at a 
residence. From tnls analysis, it is obvious 
that the impulsive character is composed of 
harmonics of tne blade passage freouency (1.16 
Hz). While the hl^iest amplitude harmonics 
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TIME. SECONDS 

Figure 1 • Time history of the sound 
pressure verletlon near 
the Hod'1 



FREQUENCY, HZ 

Figure 2 • Spectrum Analysis of 
sound variation near 
the Mod*1 


analytical method of tfc htsound coot 

Tne approach used in deveioplng the wind turbine 
sound prediction code, NTSOUND, was to apply an 
available theory used for calculating ixilse from 
conventional aircraft propellers. The 
development of such theories goes back to 1937, 
when Gutln first successfully calculated the 
noise from a propeller In a uniform flow field 
(ref. 1). Since then, Gutln' s theory has been 
exterxieu to Include the effects of non-uniform 
flow fields (ref. 2 ) and applied to helicopters 
and turbomachinery as well as propellers. The 
WTSOUNO code was written using this theory to 
provide a means of calculating sound intensity 
and freouency characteristics specifically for 


wind turblnaa in non-uniform flow fialoa, 

Tba pradofflinant sound produoad by a wind tutbina 
la aaaoclated directly with tha aarodynamic 
prasiuraa on tha bladai, Thaaa praiauroa can ba 
ralatad for convanianca to the thrust and torqua 
forcas on tha rotor. Tha thrust and torqua 
forces have componanta that are both steady and 
unsteady In time. Tha steady forcas produce 
sound called rotational noise, which consists of 
praasura variations in the acoustic field at the 
blade passing froqusncy with harmonics of 
rapidly decreasing magnitude. The unsteady 
forces may be either periodic (l.e. tower shadow 
and wlndshear) or random (l.a. gusts). Noise 
due to periodic unsteady forces may be dominant 
over rotational noise and generated hlgfier 
harmonics of amplitude comparable to that of the 
fundamental. 

The method used to determine the sound pressure 
levels in tne acoustic field is described by the 
flow Chart in Figure 3. This procedure can be 
summarized as follows: (1) calculation of the 
steady aerodynamic blade forces, (2) variation 
In these forces due to I'nsteady aerodynamics, 

3) Fourier analysis of the force variation, and 
(A) calculation of sound pressure levels in the 
acoustic field. These steps are explained In 
more detail In the following sections. 
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Flow chart of sound 
pressure level calculation 


Steady Aerodynamic Blade Forces 

The total thrust force and torque on a rotor In 
uniform flow is determined from blade element- 
momentum theory. The reader Is referred to 
reference 3 for development of this theory. 
There are several conputer programs available 
that use blade element-momentum theory. 












intJMillny the PhiTP ctiiJ# (raf. 3). The PROP coda 
raquiras mqdailnu gf tna wind turcina 
cherMctarifitlcs and oparatlng corKUtlont (i.a. 
planrernti twist, rpm, wlndopaad, ate) to 
calcuiata the staudy torque and tnruat force on 
the rotor. The tutel torque and thrust on tna 
rutur can than be rasoived into equivalent 
forces acting on each blade at a slngla point. 
The point used in this analysis is tha 73 
percent radius of the blade. That radial 
location is also used when modeling the unsteady 
aerodynsmit' variation. 

unsteady Aerodynamic (Usde Forces 

Once the steady forces have bean determined, the 
unsteady forces are obtained throu^ perturb^ 
atlon with non-cIlmensionHl force coefficients. 
Figure A srmws the the wind velocity and force 
vector diagram of an airfoil element of chord 
length c . This airfoil is operating at a pitch 
angle 0 with respect to the plane of rotation, 
The velocity of the wiftd at the rotor plane, , 
coiTblnes with the velocity due to rotation, v»i 
to give the relative velocity vector, Vr. The 
relative velocity vector acts at an angle of 
attack n with respect to the airfoil chord 
line. The relative velocity and angle of attack 
are given by 
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R 



(1> 


(3al 

■ C|^ cos(a B) + Cp sln(a + 0) ' ' 

and 

Cq • Cj^ s*n(« ♦ 0) - Cp cos(a + 6) 


Tha thrust and toraun coefficlants can now ba 
used to dsttrmine the unsteady forces associated 
with periodic variations In tha wind velocity, 

. Such e periodir variation will occur as 
Wie blade rotates through wind shear or behind 
the wind turbine support tower. The auasl> 
steady state values of blade thrust force, , 
and torque force, , at any rotor ezimuth 
position, ^ , arc given by 

arnl 

q’*' “ Q* 

c! 


(Aa) 


(Ab) 


In Which the superscript s denotes the steady 
forces and coefficients respectively. 


and 


« - - 0 


(2) 


The airfoil lift and dray coefficients are shown 
as the heavy line vnetors in Flgjre A. 



Figure 4 - Hind velocity and force 
vector dtsgram 


Fourier Analysis of the Force verlatlon 

The next step is to a perform n Fourier analysis 
of the blade force variation. The complex 
Fourier coefficients for the thrust and toruue 
forces respectively are given by 
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in wtUt^i 1 ) is the rotor speed. 

These coefficients are determined in the program 
using the IDM subroutine Fi)fUT. FORIT gives 
real Fourier coefficients Ap and Bp for the 
cosine and sine terms respectively. The 
transformation to complex coefficients is 
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( 68 ) 
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p ;> 0 


These coefficients ran be transformed .ntu a 
thrust force tioeffU' lent, Cf , acting perpendi- 
cular to top rotor i)lane and a torque force 
coefficient, Cn, acting parallel to the rotor 
plane where 


A correction to the quasi-steady state analysis 
can now t'e made by lix'ludlng the effects of 
unsteady aerodynamics. The approach used here 
to determine the response of the airfoil was 
developed by Sears (ref. A). The correction is 


VAUPATION CF W. WTSOUNO COOF 


givAD by A itlnplA AxprAtAlon CAllAd a SAArs 
runcUoo, «Mch iA miad aa a ftctar to tOA 
FourlAr coAfficiAntA. Tha ANprtAAlon uitd tu 
ApproKimAtA thA ^AfA runctlon U 
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for 


0 > 0 
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(7b) 


The coefficients An and a.p In Eduatlons (6) 
are nuitlpUea by the Sears factor before use In 
calculating sound pressure levels. 

Sound Pressure Levels In the Acoustic Field 


The mathematical relationship for calculating 
the sound pressure levels from the Fourier 
coefficients of the blade force variation mbs 
O btained by Lowson (ref. 5). The RMS pressure 
variation of the harmonic of the blade 
passage freauenuy Is given by the following 
eguatlansi 
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To verify tna accuracy of the HTSOUNP coda, 
fflaaturad data from tha Mod>l wind tuiblna wara 
uiAd, The data preaantaa hare wart taKan from 
(ieneral Flectric Company reports (unpubllAhrii tiy 
R. J. MellA. 


Correlation Mlth tna Moa>-l Haaiurad Oaia 

Oh June 10, I960 at 12:36 a.m. aouna lavala ware 
balng maaiured with a mtcropnona located about 
2A0 ft from tha md 0"1. The aina tuibine was 
operating at 34.6 rpm In a 30 npn mIdo and 
generating about 730 kw into a loan bank. Tne 
measured aouna spectrum is mown In Figure 3 ano 
has been designated case GE 160, These maasurea 
sound levels consist of botn tones and broad, 
band noise. The tone levels, mown as the 
narrow peaks, are generated by the wind 
turbine. Tne broad^iana level, mown as the 
flat areas of the spectrum, is composed of 
antiient wlno noise as well as blade vortew noise. 
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Figure 5 - Mod-1 measured sound 
levels for case G[ 180 


Since the WTSOUNO code calculates tone sound 
pressure levels, the effect of the broad-band 
noise on the spectrum itust be removed for 
comtiarlson with the code. To do this, two 
points on the spectrum were chosen that were 
believed tc be at the broad-band noise level 
(Bi and 02 in Figure 5). The broaa-band 
levels throu^out the spectrum were calculated 
using the following empirical model: 


in wnicn 

B Is the nuncer of blades 

S is the distance from tne rotor 

y« 4 are azimuth and altitude angles 
to the listener 

r^ is the blade radius where the thrust 
and toruue forces are assumed to act 

J is the standard Bessel function 

Co is tne speed of sounu 



mere . 

L' Is the broadband sound level In 
^ OB at frequency f 

L^.L- are the broa:±>ano sound levels 
' ^ in 00 it points Bi anj 

fj, f2 are the freauencles corresponding 
to points Bj end 0;^ 
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The tone-onTy levaig were then determined Dy 


*>f " 10 * - 10 ^ j (10) 


where 

is the tone aounfl level In dO et 
’ freQiiancy f 

la the measured total sound level in 
~ dB St freouency f 

Some typicHi adjusted tone*only sound presaure 
levels are shown as the sduares in Figure 5, 

Note that this correction Is small and can only 
be seen in areas where tne tone and broad>band 
levels are comparable, 

An analytical model of the Mod-1 wind turbine 
was developed for case Gt'. 180. Steady thrust 
and toroue forces, reoulred Input for the 
WTSOUNO code, were calculated using the R^OP 
performaix:c code. The wind velocity deficit in 
the wake of the tower was approximated as an 
average of the velocity profiles at the 69, 75, 
and 81 percent blacle radius. The velocity 
profiles were taken from scale model wind tunnel 
tests of the Mod-1 tower (ref. 6). Figure 6 
shows the measured wake velocities as well as 
the analytical wake model assumed. Note that 
the measured ano assumed profiles were taken at 
different distances downwind of the tower. No 
correction to the wake fur distance appears to 
be reouired, and none was made. 
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BLADE AZIMUTHAL ANGLE, DEG. 
Figure 6 • Mod-1 tower shadow 


Lowaon's eauatlon gives free-space sound 
pressures with no effect of reflection from 
nearby solid bodies, A 6 dH increase has oeen 
included in the wlbOUNtJ code analysis to account 
for reflection when using microphones near the 
ground. This correction is common practice when 
calculating airplane propeller noise at ground 
ievel (ref. 7). 

Figure 7 shows the comparison of measured and 
predicted sound levels. The squares are tne 
same as those showrt in Figure 5. The vertical 


lines are the tona levels predicted by the 
WTSOUNO code. In general tha code pradlcts the 
emplltudea of the hl{^eat harmonica vary well. 
Also tha roll-off rates (emplitude decreasing 
with fraquency) of the harmonica compare 
favorebly, The differences tnet do exist ere 
believed tb be aesocleted with the tower shedow 
model mentioned earlier. 



Figure 7 - Comparison of theoretical 
and measured sound spectra 
for case 6E 160 


To more conveniently characterlre the overall 
sound level of this spectrum, the toot sum 
sduare (RSS) of the harmonics between 20 and 50 
Hz was calculated. The lower end of this range 
was chosen as the nominal threshold of hearing. 
The upper bound was chosen because the tone 
levels above this frequency fall below the 
broad-band noise. As shown in Figure 7 tne 
20-50 HZ RSS sound pressure level given by the 
WTSOUNO code Is 2 dB lower than the measured 
data. This agreement Is felt to be very good. 

Table 1 summarizes the results from GE 180 as 
well as two other examples, cases GE 140 and GE 
lb. These twe cases were chosen to test the 
ability of the code to predict changes In rotor 
speed and distance from the machine. 

GE lAO documents the sound levels near the Mod-1 
on Aine 9, 1980 at 9:22 p.m. The wind turbine 
was operating at a reduced speed of 22.7 rpm 
xfille again generating 750 kw. The wind speed 
of 26 (tpn was lower than the 30 mph wind of GE 
160, The measured RSS sound level for GE 140 is 
17 dB less than GE 180. Note that this 
reduction is exactly the same as predicted by 
the code, 

GE IB documents the sound levels at some 
distance from the Mod-1 (3100 ft. downwind) nn 
March 31 at 12:26 pm. This location is near a 
resldefH;e that has complained of noise from the 


P? proflucing aoout 

MSO^kw into th# utility arid. Th§ miiturod 
^50 wt «ss livol la 21 dB lowar tnan oe 
IK. Tt>f prodlctofl iavti, hoMvtr» waa 27 dB 
.1^“* dctual iavai ia about 6 ee 
hi;>itr than pradlotM by nornal tpnarical 
dlaparaion, xnia 6 dB incnaaa ia btliavad to 
ba ai^latad with focualno of tna aound dua to 
^riiln and atmipnaric conditiona, Thowofi tba 
•w^t of Mpllficitlon My vary, tMa caae 
d^btratia tbo uaa of tha coot to quantify tn# 
•fftcti on aound propagation, 


Tibia 1 . Sunmary of Mod«l 
lound data cases 
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APPLICATION OF THE WTSOUND COOL 

confidence In the ability of 

*ound cbaraeterlstlcs of 
investigate 

^e Md-l problem, in particular, studies were 
^ the contribution to the sound level 
by characteristics of the flow field, and ( 2 ) 

Effect of Flow Field O^aracterlstlcs 

The STSCXe® code was used to determine the 

contribution to the souno level by the 
^aracterlstlcs of the flow throuo> the rotor 

SSInJ.rri,*?' t"* SifS?.”lSJ; 

. Shadow. As stated 

•ImiflMMi contribute 

to the fundamentii blaoe 

it esMntuii*"S^I/ "ino shear, whio 

o^ll to '^of'trlbutes 

IT ^ i?. . •'•hWhics. However witn 

Shadow, mny hamonics an 
•udlbl^fSJ^^^"'* shplltudes up into fS 

c^trlbuti^^I^L”'^?* Predomlnent 

problem of Mod-1 is 
tower*"'* deficit In the wake of the 



FREQUENCY, H2 

Figure 8 - Effect of wind shear and 

tower shadow on wind turbine 
sound level (Mod-1 case GC IS) 

Effect of Operating Conditions 


unflar i rh^K. cne conditions 

j^er «e>icn the highest sound levels ere 

effe^R^^Af^''? JITSOUNO code was used to study tne 
effects of various operating parameters, in 
effects of power, rotor speed 
wind speed, and direction were lAvestigstSS ’ 

ngure 9 snows the predicted sound levels versus 

SrjoMnI°L,S“‘" 23 rpm operatio" ' 

oparatlon as tne 
fP*M «rles from cutin through rated to 
cutout wind speed (6, is, and 20 m/sec 

Analysis shows about a 12 

speeds for both 35 and 23 rpm operation. 
^Is agrees reasonably well with enplrlcal 
modeling of the measured data that Indicates a 
variation of about la oe over those wind speeds. 

SSlS 23 rpm, the sound 

levels predicted In Figure 9 are about 11 dB 

Sit'.Kl!'’.!?'" ">"■ .q»» 

"•■sured data taken in June 1980. 
tkjring those tests, the Mod-1 was ooeratej at 

Sse“datS*2hoJn‘“ • 

rsfo!*2e.’5‘‘o"?'*2; 


SSiSSa^?® ih the wiis foT:?™'*" 

OpAfaiion Ih this 

region Is chirecterlreo ai off-optlmum under 
^Is Situation the pitch of the bleSTis 

SS^JtiSo^Tl^ Taatner, to limit power below 
ihLJ machine. The shaded ares 

levels fflty be produced 

to be suOMrte^h conclusion seems 

laoT75Q*k»*^*^5r^\^“'* ogeretion than Cf 

ISO 1750 kw, 26 itiph) end Is hijy^er In souno 
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than tha 11 (£ dlffarnnce assoclateo with 
rpm only, 


Th« iouna levaU at off-optKnuin oparatlon abova 
ratad hava not baan datarmlnad. Iinprovainanta In 
tha unataady atrodynamlc forca calculation in 
tna ttTSOUNd code ara balng incorporatad to 
anatyia thia arta of oparatlon, 
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Figure 9 • Effect of power, rpn, and 
windspeed on Mod-1 wind 
turbine sound levels 


The wTSOUND code can also be used to calculate 
the olrectivity pattern of tne sound. Figure 10 
shows the directivity pattern In a plane 1100 
ft. below the Mod-1 hub height. The direction 
of highest sound pressure is directly downwind; 
nowever, it Is nearly as ni^ upwind. Minimum 
sound generation occurs in the rotor plane where 
levels are about 16 dB lower than downwind. 

This compares reasonably well with measured date 
in the near field showing about 16 dB variation. 


OONCLUSIONS 

The WTSOtiND computer code shows generally good 
agreement with sound spectra measured in the 
vicinity of a wind turbine. In the far field, 
however, correlation of the absolute amolltude 
of the sound level is complicated by propagation 
effects. For the case in this study, terrain 
and meterologice 1 conditions caused an increase 
of about 6 db. 

Analysis using the SOUND code shows that the 
predominant contributor to the noise problem of 
Mod-1 is the wind velocity deficit in the wane 
of the tower. Changes in the aerodynamic 
forces, as the blades pass through the deficit, 
produce sound pressure variations in the 
acoustic field. 

The level of the sound pressure variations are 
most directly affected by rotor speed and 
windspeed. Reducing the rotor speed from 36 to 
23 ipm is predicted to reduce sound levels by 


■bout 11 dB, The increase In sound levels with 
windspeed is predicted to oe 12 dB between cutin 
end reted, 
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Figure 10 - Predicted sound directivity 

S Bttern In a plane 1100 ft. 
elow the Mod-1 
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